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Abstract 

We investigate new and unusual signals that arise in theories where dark matter 
is asymmetric and carries a net antibaryon number, as may occur when the dark 
matter abundance is linked to the baryon abundance. Antibaryonic dark matter can 
cause induced nucleon decay by annihilating visible baryons through inelastic scattering. 
These processes lead to an effective nucleon lifetime of 10 29 — 10 32 years in terrestrial 
nucleon decay experiments, if baryon number transfer between visible and dark sectors 
arises through new physics at the weak scale. The possibility of induced nucleon decay 
motivates a novel approach for direct detection of cosmic dark matter in nucleon decay 
experiments. Monojet searches (and related signatures) at hadron colliders also provide 
a complementary probe of weak-scale dark-matter-induced baryon number violation. 
Finally, we discuss the effects of baryon-destroying dark matter on stellar systems and 
show that it can be consistent with existing observations. 



1 Introduction 



Cosmo logical observations indicate that about 4.6% of the energy density of the Universe 
consists of baryonic matter, while 23% is dark matter (DM) [TJ. Neither of these results 
can be explained with our current understanding of elementary particles, the standard 
model (SM). Cosmology therefore requires new fundamental physics, and it is important 
to find ways to detect such new physics experimentally. 

In the majority of new-physics scenarios, the generation of baryons and DM occurs 
through unrelated mechanisms, offering no explanation for the similar magnitudes of their 
cosmological densities. The most thoroughly studied scenarios involve baryon production 
from CP-violating non-equilibrium processes during the electroweak phase transition, from 
decays of right-handed neutrinos, or from the coherent evolution of scalar fields [2] , while the 
DM relic density is determined by thermal freeze out when a non-relativistic stable species 
falls out of equilibrium [3]. In this context, there is no reason to expect similar cosmological 
densities of baryons and DM. 

This apparent coincidence may instead be a clue that both types of matter have a common 
origin. Several models of asymmetric dark matter (ADM) have been proposed along these 
lines where the DM density carries a net (approximately) conserved global charge shared by 
theSM |H El El El [H El HH HH H31 [Til [15], such as baryon number B. These models 
generally fall into two classes depending on how the charge asymmetry is created: 

1. An initial charge asymmetry, generated in either the visible or DM sector, is partitioned 
between the two sectors by chemical equilibration through a transfer operator [U [9]. 
These charges are separately "frozen in" once the transfer operator goes out of equi- 
librium. 

2. Non-equilibrium dynamics generate equal and opposite charge asymmetries in the vis- 
ible and DM sectors, without any net overall charge asymmetry pTJl [TTl [12 l [T3l [T^l [15]. 
In order to avoid washout, transfer operators must always be out of equilibrium once 
the asymmetries are created. We term this process hylogenesis ("hyle" = matter) [13] . 

Subsequently, in both cases, DM particles and antiparticles are assumed to annihilate effi- 
ciently, leaving only a remnant asymmetric component determined by the charge density. 

In the present work, we investigate novel experimental signatures from hylogenesis scenar- 
ios where DM carries B [13]. Here, the lowest- dimensional, gauge- invariant transfer operator 
is given by 

^-Iu^rHIix. (1) 

where k label generation. In this case, DM has two components, a fermion/scalar pair 
(\l/,$) with total baryon number Bq + B$, = — 1. Stability of both ^ and $ requires 
{my — m$| < (trip + m e ). This operator is inactive cosmologically and does not wash out the 
baryon-DM asymmetry for A > 100 GeV provided the asymmetry is created at relatively 
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low temperatures, below about a GeV. Our results may also be applicable to other ADM 
scenarios in which the same transfer operator appears with one or both of $) making up 
the DM. 

In hylogenesis, the Universe is net S-symmetric, and therefore the baryon asymmetry 
carried by DM is equal and opposite to that in visible baryons. Specifically, in the hylogensis 
scenario of Ref. [15J, this has two important consequences: 

• The number densities of $) satisfy n# = n<$> =Ub- Therefore, cosmological ob- 
servations imply (m^ + m^)/m p = Q DM /Q b rs 5. Together with the DM stability 
requirement, we have $ « 1.7 — 2.9 GeV. 

• There exist many scenarios in which ADM is coupled to the SM via the "neutron portal" 
operator undndn pE], generally falling into the first class of "chemical equilibration" 
scenarios. In this case, DM is very often baryonic (Bdm > 0), while for hylogenesis 
scenarios the DM must be antibaryonic (Bdm < 0)lll 

In Ref. [15] , we presented a specific realization for hylogenesis where the operator in Eq. (TjQ) 
arises by integrating out heavy Dirac fermion mediatiors Xi t 2- Out-of-equilibrium decays 
of X\ during reheating generate equal-and-opposite dark and visible baryon asymmetries. 
Furthermore, in this scenario $) are charged under an additional hidden U(l)' gauge 
symmetry, that couples to the SM via kinetic mixing with hypercharge, to facilitate anni- 
hilation of the symmetric DM densities. More details about this realization are given in 
Appendix [A] 

Hylogenesis models with the operator of Eq. (CQ) have an interesting and unique signature: 
antibaryonic DM particles can annihilate visible baryonic matter, termed induced nucleon 
decay (IND). IND is a novel and unusual DM signal. These events are — > and 
— > \I/M, where N = (n,p) is a nucleon and M is a meson. Since the DM states are 
invisible, this mimics nucleon decay, with an effective lifetime dependent on the local DM 
density. This signal offers the new and exciting possibility of searching for DM in nucleon 
decay searches in deep underground detectors such as SuperKamiokande [IS] . 

In Section [21 we compute the rates and kinematics of IND, and discuss the implications 
of IND for existing nucleon decay searches. Our main conclusions are: (i) due to different 
kinematics of IND, compared to standard nucleon decay, existing bounds do not apply over 
most of the region of parameter space of our DM model, and (ii) the effective nucleon lifetime 
can be around 10 29 — 10 32 years, if the new physics scale A in Eq. ([I]) is the weak scale. 

Hadron colliders can probe (anti)baryonic DM scenarios through direct production of the 
weak-scale mediators transfering B between visible and dark sectors, discussed in Section |3j 
The operator of Eq. (CQ) can give rise to observable signatures in the form of monojets (or 
jets and missing energy). Within our hylogenesis model, if IND were detected in nucleon 
decay searches, monojet signals are inescapable and should be observed at the Large Hadron 
Collider (LHC). 

1 In supersymmetric models, light superpartners can play an important role in chemical equilibration [16] . 
potentially affecting the sign of -B^m m equilibration ADM scenarios. Also, for other applications of the 
neutron portal operator to baryogenesis, see Ref. [17] . 
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Figure 1 : Feynman diagrams for induced nucleon decay. Box denotes IND vertex from Eq. ^ . 
Circle denotes strong interaction vertex given by Jz?o in Ref. [23]. 



IND processes can potentially be relevant in stars. Since the IND rate scales with the DM 
density, capture and accretion of DM in stars is important, possibly leading to modification 
of stellar evolution due to baryon destruction and energy injection. In Section HI we consider 
these effects in neutron stars, white dwarfs, and main-sequence stars. Our conclusions are 
summarized in Section [H 



2 Dark Matter Detection in Nucleon Decay Searches 

Dark matter particles $) can annihilate nucleons N, producing an energetic meson M 
through inelastic scattering 

$iV — > #M , — ► $ f M , (2) 

shown in Figure [U We restrict our attention to single meson final states, although multi- 
meson events are allowed and may have comparable rates. In general, both down-scattering 
and up-scattering can occur, defined as whether the heavier or lighter DM particle, respec- 
tively, is in the initial state; however, the latter is kinematically forbidden if |m# — m$| > 
mjy — itim- Since neither the initial DM nor final anti-DM particles are observed directly, 
these processes mimic standard nucleon decay events N — > Mv [19], with an undetected 
final state neutrino v (or antineutrino v). 

The observable meson energy for each case in Eq. (j2J) is, respectively, 

(m N + m$) 2 + m M - m% (m N + m*) 2 + m 2 M - m% 

&m — —, ; ; , &m — 7T{ ; \ > i 6 

2{rriN + m<i>) 2{mN + m^) 

neglecting the initial kinetic energy of the DM particles (vdm ~ 10 -3 ). 

Although hylogenesis (as a baryogenesis mechanism) works for any quark flavor, the 
most interesting signatures arise if the IND operators involve u, d, s quarks only. Restricting 
our attention to the lightest mesons (ir, K, 77), the final state meson can have much more 
kinetic energy in IND than for standard nucleon decay (SND), summarized in Table. [TJ For 
fixed masses and m$, the meson is either monochromatic or bichromatic (up to Fermi 
motion), depending on whether up-scattering is forbidden or allowed. The range of momenta 
Pm corresponds to the allowed range m^^ ^2 — 3 GeV, provided that $) are stable and 
account for the observed Q,dm- 
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Decay mode 


„SND 

Pm 


< D [up] 


Pm° [down] 


t% nd bound (xlO 32 yr) 


N -> 71 

N ^ K 
N -> 77 


460 
340 
310 


< 800 

< 680 

< 650 


800 - 1400 
680 - 1360 
650 - 1340 


t* nd > 0.16 [2D] , t$ nd > 1.12 [21] 
rf™ > 23 [IB] , Tf™ > 1.3 [IS] 
rf^ > 1.58 [21] 



Table 1: Comparison of meson M = (tt,K,7]) momentum pm (MeV) for standard nucleon decay 
(SND) and induced nucleon decay (IND) from DM, for up- and down- scattering. 

2.1 Nucleon decay searches 

Existing searches have been optimized for meson momenta pff D ~ 300 — 450 MeV, while 
for IND mesons are typically much more energetic, with momenta plf ~ 1 GeV. Here, we 
briefly summarize existing SND search strategies and how they might be adapted for IND 
searches. 

p — > K + u, n — > K°v: The Super-Kamiokande experiment, a water Cerenkov detector, 
provides the strongest limits on these channels. For K + , they have three searches: (i) 
K + — > 7T + 7r°, giving three Cerenkov rings, (ii) K + — > fi + with a prompt 7 (from 16 O — > 
15 N* — > 15 N + 7), and (in) a mono-energetic /i + from K + — > with no prompt 7. All 
three searches assume, as is the case for p^+ ~ 340 in SND, that the K + is emitted below 
Cerenkov threshold (f3 < 0.75) and comes to rest before decaying. For IND, we estimate that 
an 0(1) fraction of i^ + 's do come to rest before decaying. However, except for up-scattering 
events close to kinematic threshold, the K + from IND has (3 > 0.75, adding an extra ring to 
the event topology. Furthermore, this additional radiation may make finding the prompt 7 
in search method (ii) more difficult. For K°, they have two searches: (iv) K° s — > 7r°7r° — > 47, 
giving four e-like rings, and (v) K° s — > 7r + 7r~, giving two /i-like rings. Both searches assume 
200 < pxo < 500 MeV, thereby excluding IND events (again, except for up-scattering near 
kinematic threshold). One difficulty in search (iv) is identifying all four e-like rings. For 
IND, this may be more difficult as the rings would be more overlapping due to relativistic 
beaming. On the other hand, search (v) is promising for IND. In SND this mode suffers from 
a small efficiency that the 7r are both above Cerenkov threshold. With greater energetics 
in IND, the efficiency may be much larger. 

p — > 7i + v: The best limit is provided by the Soudan 2 experiment, an iron tracking calorime- 
ter. Nucleon decay event candidates were required to have a single 7r + track, with ionization 
consistent with mass or m M , initial momentum 140 < p 7T + < 420 MeV, and visible 
endpoint decays (tt + — > fi + — > e + ). Their simulations showed that a tt + originating from 
within an iron nucleus loses on average half its initial momentum. At higher p n + , IND events 
may be more visible due to reduced background from atmospheric neutrinos. However, it is 
unknown to us what is the average momentum deposition in iron of the tt + at much higher 
energy, and whether this can lead to fragmentation of the parent nucleus into exotic nuclear 
states. 

n — > 7r°i>, n — > f]u: The best limits on these modes come from the IMB-3 experiment, a 
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water Cerenkov detectorH The 7r° — > 77 channel may be more difficult at higher energies: 
due to decreased separation angle of the two photons (from relativistic beaming), they can 
appear as a single electron-like track. The n — > rju, 77 — > 77 channel will have greater photon 
separation and may be more promising. We find that the IND rates into tt° and 77 final states 
are comparable and are sensitive to the same underlying IND operator (shown below). 



2.2 Effective nucleon lifetime from IND 

An effective IND lifetime can be defined as the inverse scattering rate per target nucleon, 
r^ 1 = n DM (av) IND , with local DM number density n DM = pDAi/( m ^ + m <s>) and IND 
scattering cross section {av) IND . Numerically, we have 

« (10- yrs)"' x ( n J°» ) ( . (4) 



0.3GeV/cm 3 y V10~ 39 cm 3 /s 

Next, we compute (ctv)ind using chiral perturbation theory. We perform an expansion in 
powers of Pm/(47t/), where / ~ 139 MeV is the pion decay constant, and truncate at leading 
order. Since for IND we expect pm ~ 47r/ ~ 1 GeV, our calculations should be regarded as 
order-of-magnitude estimates at best. Our analysis closely follows SND rate computations 
in Ref . 



There are four effective interactions that are relevant for IND processes with single meson 
final states. These are given by Jzf int = J^CjOj, with operators (given in two-component 
spinor notation) 



O^e^rt^s) (5) 

2 = ^ e^*[(ds4)(«£*«) + (*£«£) (<W - 2«4)(4^)] (6) 

3 = ^e^*[(dS4)(«M - (7) 

where 7 are color indeces, and the coefficients q have mass dimension —30 The linear 
combinations have been chosen such that have strong isospin I = (|, 0, 1), respectively. 
Here, it is useful to write Jz^ nt = Tr(cO) where 



eg. _ ca_ Q 

'2 



ci -\/^ C2 / 



2 We note that the IMB-3 experiment found an excess in events with total energy 900 — 1100 MeV (20 
events vs. 6.1 expected background) [22] ■ Nearly all these events had between 2 — 4 Cerenkov rings, large 
missing momenta (400 — 1100 MeV), and large invariant masses (600 — 1100 MeV). In IND, a large missing 
momentum would be expected, while a large invariant mass could arise through heavy meson or multi-meson 
final states. 

3 A fourth operator e a p-y^(sRUji)(sii^f u) is relevant only for multi-kaon final states. A fifth operator 
e a/ 3 7 md R s^)(u 7 R ^ R ) + (s%u p R )((P R * R ) + (u R d R ){s R $ R )} vanishes by a Fierz identity. 
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Figure 2: Induced nucleon decay cross sections (o~v)ind for p,n — > 7r + ,7r° (left) and n — >■ r] (right) 
as a function of fermion DM mass mq, for \ci\ = TeV~ 3 . Dotted (dashed) lines denote — > ^fM 
(N 1 ® -> &M). Solid lines denote total rates iV<£> -)■ + -)■ $+M. (av) IND = 1CT 39 cm 3 /s 
corresponds to lifetime t^ nd = 10 32 years. 



with = (u,d,s)ft. Under SU(3)l x Sf/^R chiral symmetry transformations, the right- 
handed quark fields transform in the (1, 3) representation, qn — >■ -Rg/j (where i? G ^[/(S)/?), 
while the IND operator transforms in the (1, 8) representation, O — >■ ROR) . If we treat 
c as a spurion in the (1,8) representation, Jzfi nt is invariant under chiral transformations. 
(The DM fields (\1/, $) are chiral singlets.) 

The IND interactions of baryon and (pseudo-Goldstone) meson fields are determined by 
the chiral transformation properties of the spurion c. Following the conventions of Ref . [23J , 
the only invariant operator is 



Sf IND = fiTr[c^{B R ^f R )^] 
where £ = exp(iM/ f). The meson and baryon fields are 



(9) 



/ JL + *£. 



M 



71 



V 



7T 



\ 



K c 



B 



AH 

V6 



P 



• (10) 



3 JV / 



Eq. (jUJ) is invariant since the quantity ($Br£) is in the (1,8) representation [23]. The 
unknown overall coefficient (3 = 0.014(1) GeV 3 has been computed using lattice methods 



The Feynman diagrams for IND are given in Fig. [IJ It is straight-forward to derive the 
Feynman rules for the interactions of baryons, mesons, and DM by expanding the matrix 
expressions in Eq. (Q and working to linear order in 1/f. From these, we compute the 
matrix elements and cross sections for IND processes. We assume that only one coefficient q 
is non-zero at a time. N — > it, i] modes depend only on ci, while N — )■ K modes are governed 
by c 2 , 3 - 

Our numerical results are shown in Figs. [2] and [3j for the five decay modes we consider. 
We plot the velocity-weighted cross sections (av)jND for each mode, as a function of ss 
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Figure 3: Induced nucleon decay cross sections (o~v)ind for p — > K + (left) and n — > K° (right) as 

a function of fermion DM mass mq, for | C2,3 1 = TeV~ 5 . Dotted (dashed) lines denote — > $K 

(Nfy — > &K) from operators 02,3- Solid lines denote total rates — > ^>K + — > Grey 
regions show where existing nucleon decay bounds apply, described in text. 



(5m p — m$), for the allowed range 2m p < mq, < 3m p (as required by DM stability). In 
Fig. El the solid lines show the total cross sections for the channels p — > tt + and n — > tc° (left) 
and n — > r\ (right). All three modes arise from the same IND operator Ox, and we have fixed 
ci = TeV~ 3 . The individual cross sections for — > (iV$ — > M^f) are shown by the 

dotted (dashed) curves. Although the individual rates vanish where kinematically forbidden, 
the total rate is always non- vanishing. For — m$| < rriN — itlm, both up- and down- 
scattering rates are non-zero, and IND is bichromatic. Otherwise, only down-scattering is 
allowed, and IND is monochromatic. 

In Fig. [HI we show the strange IND channels p — > K + (left) and n — > K° (right). Here, 
there are two relevant operators 02,3. The solid lines denote the total IND cross sections 
induced by each Oi independently, taking q = TeV -3 , for i = 2,3. The individual cross 
sections for N^/ — > (N& — > M\]/) are shown by the dotted (dashed) curves. Again, 

the kaons can be monochromatic or bichromatic, depending on whether up-scattering is 
kinematically allowed. 

Existing nucleon decay bounds do apply in select regions of parameter space for up- 
scattering close to threshold, where the meson momentum is reduced. We illustrate these 
regions in Fig. [HI shown in grey, for the case of Super-Kamiokande. For p — > K + , this region 
corresponds to (3k+ < 0.75: in this case, the K + is below Cerenkov threshold and the event 
topology is identical to p — > K + v. For n — > K°, this region corresponds to a kinematic 
window 200 < px° < 500 MeV, as in their n — > K°v search. We emphasize that Super- 
Kamiokande bounds constrain only the up-scattering IND rate, which can be suppressed 
compared to down-scattering in these parameter regions. 

The total IND rate depends sensitively on the unknown mass scale Aj^r, = |q| _1//3 . The 
total IND cross sections and nucleon lifetimes for all channels are comparable, scaling as 

/ \ m-39 3 / ( A-INd\ 6 -, n 32 ( ^IND \ 6 ( PDM \ 

(av), ND « 10 cm /s x j , r N M 10 vr x j ( ^jGeV]^ ) (U) 
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As we show in Sec. [31 the collider bound on this scale is Ajnd > 300 GeV. Therefore, tn 
can in principle be as low as 10 29 years. It is likely that such a short lifetime would be 
excluded from nucleon decay searches, but no dedicated IND search has yet been performed. 
We also note that future nucleon decay experiments are envisioned to have markedly better 
reach for rjv- For example, a water Cerenkov detector with 10 4 kton x year of exposure can 
reach r/v ~ 

10(34-35) 

yr [251 12S] ; a similar capability is expected to be achieved with a liquid 
Argon detector with 10 3 kton x year of exposure [27]. For kaon final states, liquid Argon 
technology is expected to provide improved efficiency due to better imaging capabilities [27J. 

Lastly, we note that the IND cross sections satisfy certain relations, as a consequence of 
strong isospin symmetry: 

(avf^o = 2 (ov)%£ , M&£ + = (H?^f • (12) 
The latter relation holds only if iV — > K modes are dominated by either O2 or O3, as assumed 
in Fig. [31 If neither operator is negligible, then (cu)?^ 7^ ( au )?iviT '■ m this case, both 
K modes are complementary and can be used to disentangle the underlying IND operator 
structure. 



3 Collider Signals from Hylogenesis 

ADM scenarios rely on transfer operators to connect global charge between the visible and 
dark sectors, and these operators can be probed at high-energy colliders. In the specific 
hylogenesis model of Ref. [15], B is mediated between the two sectors by heavy Dirac fermions 
Ai,2 (ttixi < m x 2 )i through interactions of the form 

- * D E % (44)(<z4) + (a (X a , L * L + X a , R * R )* + h.C. (13) 
a=l,2 

where i, j, k label generation, color indices are implicitly contracted antisymmetrically, and 
other fermion contractions are also possible. Integrating out X 12 generates operators of the 
form of Eq. (JTJ) . 

The X12 particles can be produced at high-energy hadron colliders through the operator 
of Eq. (pp. With decays X l>2 this gives rise to events involving missing energy and 

one or more jets. In this section, we investigate the sensitivity of the Tevatron and the LHC 
to such events and we derive a corresponding bound on the heavy mass scale M suppressing 
the neutron portal operator. Related studies in the context of WIMP and other dark matter 
candidates can be found in Refs. [281 1221 US EH E2l [33] . 

To be concrete, we will focus on the lighter state X\ = X, with the specific interaction 

-se d ^(xls R )(u R d R ) + cxm + h.c, (14) 

We expect other flavour structures and fermion contractions to give qualitatively similar 
results. The operator of Eq. (I14p can give rise to processes of the form 

q{pi)q'{p2) -> fipa) *(p 4 ) $ f (p 5 ) (15) 
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through either a real or off-shell X, where q, q', q" = u,d,s quarks. The corresponding 
summed and averaged squared matrix element can take two possible forms, depending on 
how the fermions are contracted. They are: 

2 

(Pi • P2) [2(p3-?)(P4-g) - (g 2 - rn 2 x )(pi'P4)\ ; s-like 
\M\ 2 ={ (16) 





AC 
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1 




M 2 ' 




q 2 —m%+iT x m x 


ll 


AC 


2 


"la; 




I? 2 ' 




q 2 ~m x +ir x m x 



2 



(pi -Ps) [2(p2-g)(p4-g) - (g - w4-)(p2-P4)] ; Mike 

Here, g = (P4+P5) = (P1+P2— P3) is the momentum carried by the intermediate X state, and 
T x = £ 2 mx/167r is the width of the X state which we assume decays mainly into The 
s-like form correponds to the case where both initial state fermions are contracted together 
in the underlying operator while the t-like form corresponds to a contraction between initial 
and final state quarks. 

At the Tevatron and the LHC, we find that the full cross section derived from the 
matrix elements of Eq. (Tl6|) is frequently dominated by the pole in the intermediate X 
propagator. This corresponds to the production of an on-shell X state whose mass is not 
much smaller than the higher- dimensional operator scale M, and corresponds to a large 
momentum transferal In this case the higher-dimensional operator structures we are using 
become unreliable and the full dynamics of the unknown ultraviolet completion becomes 
relevant. Without specifying the underlying theory, we can still parametrize the generic 
behaviour in a reasonable way as follows. For s-like contractions, whose structure would 
arise most naturally from a boson in the s-channel, we make the replacement 

— A 

where s is the parton-level Mandelstam variable and T is the decay width of the mediator. In 
the absence of an underlying theory, we parametrize this quantity as T = CM and consider 
the values C = 1/5, 1/50. For t-like contractions, we make the replacement 

W T=Mi> (18) 

where t = (pi — P3) 2 is the parton-level Mandelstam variable. Since t is negative, we do not 
bother adding a width term. 

We apply these matrix elements to compute the leading-order (LO) monojet production 
cross sections at the Tevatron by convolving with CTEQ6.1M parton distribution func- 
tions [31] and integrating over phase space. To match the most stringent Tevatron monojet 
search bounds, we impose a cut of pr > 80 GeV and < 1.0 on the outgoing jet. Following 
Ref. [30], we also apply a flat efficiency factor of 40% to connect our parton-level cross 
section to the full hadronic jet reconstruction at the Tevatron detectors. The cross sections 
computed in this way (after applying cuts and the efficiency factor) are shown in Fig. H] for 
the parameter values mx = 0.75M, 1.5M and V = M/5, M/50. We also set A = 1 and 



This is a necessary condition for hylogencsis to create a sufficiently large asymmetry. 
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Figure 4: Leading-order monojet production cross sections at the Tevatron subject to the 
cuts described in the text. We show lines for mx = 0.75M, 1.50 M and T = M/5, M/50, 
and we set A = 1 and ( = 0.7. 

( = 0.7. This figure shows a significant resonant enhancement when the intermediate state 
is narrow. A similar enhancement was seen in Ref. [32J. The resonant enhancement only 
becomes fully operational at M > 300 GeV due to the cut on jet pt- The horizontal dashed 
line in Fig. H] represents the current Tevatron 2a exclusion limit on the net monojet cross 
section of 664 fb (after cuts) [30], based on the CDF analysis of Refs. [35j [36] which uses 
the same set of jet cuts as applied to our signal estimates. This limit translates into a lower 
bound of M — 200 — 700 GeV for the operator of Eq. ( fT4l) . depending on the mass of the X 
state and the width of the unspecified intermediate state. 

Monojet signals can also be detected at the LHC Due to the expectation of a significant 
amount of associated QCD radiation, Ref. [37] investigated the reach of an inclusive search 
for a hard jet plus missing energy search at ATLAS with no veto on additional hard jets. 
To match this analysis, we compute the inclusive leading-order parton-level cross section 
at the LHC with y/s = 14 TeV subject to the cuts pr > 500 GeV and \r]\ < 3.2 on the 
outgoing jet. We also rescale the cross section by a conservative acceptance/efficiency factor 
of 85% [30J. The corresponding cross sections are shown in Fig. [5] for several values of M 
with m x = 0.75M, 1.5M and T = M/5, M/50. We also set A = 1 and ( = 0.7. Based on 
the background analysis of Ref. [37], the study in Ref. [31)] estimated a net SM background 
production rate after the applied cuts of &bg — 200 fb. Applying a simple S/ \/l3 > 5 
measure on the detection significance, this leads to a sensitivity to monojet cross sections as 
small as 70 fb (7 fb) with 1 fb' 1 (100 fb' 1 ) of data at = 14 TeV. This is shown by the 
dotted lines in Fig. [5j The resulting LHC reach for the operator of Eq. ({14"]) lies in the range 
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Figure 5: Jet plus missing energy production cross sections at the LHC (14 TeV) subject to 
the cuts described in the text. We show lines for nix — 0.75 M, 1.50 M and T = M/5, M/50, 
and we set A = 1 and ( = 0.7. 

M = 1 — 4 TeV. This is competitive with the sensitivity to this operator from searches for 
IND in existing nucleon decay experiments. 

Additional quark operators beyond that given in Eq. ( 1T41) and considered above are 
expected to yield qualitatively similar collider rates and signatures in most cases. An 
interesting further possibility are quark operators involving top or bottom quarks. With 
a bottom quark in the final state, the monojet signal could be augmented with a fe-tag. In 
the case of a top quark in the final state, the signal would be a hard single top quark with 
large missing energy. While the search reach in monobottom and monotop channels is likely 
greater than for the light quark channels, the operators involving third-generation quarks do 
not correlate directly with IND processes. 

The hylogenesis model in Ref. [15] also contains a hidden U(l)' gauge symmetry that 
couples to is spontaneously broken at the GeV-scale, and couples to the SM only 

through kinetic mixing with hypercharge. Direct pair production of DM particles via the 
corresponding Z' vector boson with an associated jet can be another source of monojet 
signals [301 El]- However, since this vector is relatively light, the existing bounds from the 
Tevatron and the expected reach of the LHC are both much weaker than the limits from 
DM direct detection through elastic scattering mediated by the Z'. 
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4 IND and Stellar Evolution 



DM can be captured in stars by elastic scattering with the nuclei they contain. Once cap- 
tured, ordinary self-annihilating DM will collect in the middle of the star and annihilate with 
other DM particles there, releasing energy [381 EHJ HO] ■ Self-annihilation is not possible for 
hylogenic dark matter (hDM) due to its conserved global charge. Instead, it can potentially 
destroy baryons within the star via IND processes such as ^ N — > If it remains within 

the star, the anti-hDM $t reaction product can then annihilate with a $ particle already 
captured in the stellar core. The net result of this chain is the destruction of a baryon and 
the release of energy, both from the meson decay and the annihilation step. 

In this section we investigate the effects of IND processes on several varieties of stellar 
species including neutron stars, white dwarfs, and main-sequence stars. To be concrete, we 
will concentrate on the specific model of hylogenesis presented in Ref. [15] and described 
in Appendix |X] with IND mediated by the operator of Eq. ([1]). We will assume a fiducial 
spin-independent proton scattering cross section of a^ 1 = 10 -39 cm 2 (and erf 1 = 0) for both 
\1/ and $ induced by their coupling to the kinetically-mixed U(l)' vector bosonjf] as well as 
an hDM-anti-hDM annihilation cross section of (crv) ann = 10 _25 cm 3 /s. We shall consider 
two cases for the IND cross section: a large value of ((Jv)in£) = 10 _39 cm 3 /s, and a small 
value of (ctv)jnd — 0. The general behaviour for other IND rates will lie somewhere between 
these two extremes. 

4.1 Stellar Capture and Annihilation 

Relic \1/ and $ particles will be captured in stars by scattering with nuclei to energies below 
the local escape velocity. Once captured, an hDM particle will undergo further scatterings, 
thermalize with the baryons in the star, and collect within the stellar core. This occurs 
quickly relative to the lifetimes of the stars we consider here for our large fiducial value 
of the proton scattering cross section. Once they thermalize, the DM particles are largely 
confined to a core of radius [391 HOI HI] 



where i = ^, $, T c is the mean temperature and p c is the mean (baryon) density in the 
stellar core. 

The evolution of the total numbers of $, and $^ within a star is described the 
5 For DM masses below 3 GeV, this value is consistent with existing direct detection searches p!TlFr21H5] . 




(19) 



13 



following system of equations: 



dt 
dNg, 

dt 
dN<s, 

~dT 
dt 



- AyN 9 N$ - B^Ny (20) 

- AyNyNy + e^B^N® (21) 
C$ - A$N N& - B^N* (22) 

- A$N*N& + e&BqN* (23) 



Here, the Cj coefficients are the hDM capture rates, the A4 coefficients describe hDM-anti- 
hDM annihilation, and the Bi coefficients describe IND. A general expression for Cj is given 
in Ref . [40] , while the A and B coefficients are given to a good approximation by 

Ai ~ (ov\ ann j (47^/3) , (24) 

Bi ~ (crv)u ND {p c /m n ), (25) 

where m n is the mass of a nucleon. The terms appearing in Eqs. (12"T1) and fl2"3"|) are the 
probabilities for the anti-hDM products of IND to be captured by the host star after they 
are created. In certain regimes additional processes can influence the evolution of the stellar 
populations of hDM and anti-hDM such as evaporation [39, 45J and direct annihilation to 
baryons (e.g. *S> $ — > N M). We will discuss these effects when they may be relevant. 



4.2 Neutron Stars 

Neutron stars are very dense objects supported by the Fermi degeneracy pressure of their 
neutrons. Despite their name, they also contain a significant mass fraction of protons and 
heavier nuclei, and the nuclear state of their cores is not fully understood |46j . Typical 
neutron star parameters are mass M = 1.4M , with M ~ 2.0 x 10 30 kg the solar mass, 
radius R = 10 km, core temperature T c = 10 5 K, and core baryon density p c = 1.4 x 
10 18 kg/m 3 [3H HE nuns]. 

The rate of capture of ^ or $ upon a typical neutron star (including general relativistic 
corrections) is [471 14H1 149] 

C, c 2.5 x 10*V ( 1 ^L\ (^-) (™EE!l) f, (26) 

where pdm is the local DM energy density, v is the local DM velocity dispersion, and 

/ = m in {l, (x p a p + x n a n )/(2 x 10~ 45 cm 2 )} , (27) 

with x p and x n being the proton and neutron mass fractions. The factor / accounts for the 
saturation of the cross section at the cross-sectional area of the star. This saturation sets 
in when the star becomes optically thick to DM - when a DM particle impinging upon the 
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star is likely to scatter multiple times with nucleons in the star. For the fiducial nucleon 
scattering cross section we are using and assuming a proton mass fraction of x p = 0.1, we 
find that neutron stars are optically thick to both \1/ and $. This leads to Cq, = C$ = C as 
well as = 1 = e$. We also find 

,fT c \ 1/2 /3 GeV\ 1/2 / 1.4 x 10 18 kg/m 3 
r,,*~ (140 cm) — —— . » ( 28 ) 



io 5 Ky V m i 7 

implying A* ~ 5 x 10 _32 s _1 and Bi ~ 0.9 s" 1 (0 s -1 ) for the large (small) IND rate (ctv)jnd 
10~ 39 cm 3 /s (0cm 3 /s). 



Consider first the case of a large IND rate. Using the evolution equations of Eqs. 
and assuming negligible initial hDM densities, we find that the numbers of hDM and anti- 
hDM particles within a typical neutron star reach a steady state. This behaviour is shown 
in Fig. Ofor the fiducial cross sections given above and = 2.85 GeV, m$ = 2.05 GeV. 
The steady-state populations in the limit of Bi 3> A4 are approximated well by the analytic 
expressions 

C AT B * (A±l_\ AT g AT B * (0^ 




*-£?*(i+0' w ~ a, vi + rv ^(1+e- 1 )' *' v i+e 

with £ = AyB$/Ag,B 9 E 

A steady state is also attained when one of the IND rates vanishes 
due to kinematic suppression. For example, with B^ 3> A^ ; $ and B$> — > 0, we find 

iV » = ^ ** = °> ^=±^- + 1/(^-1 +-f". ( 3 °) 

where in the last equality the plus sign corresponds to $ and the minus sign to $ f . The 
solution for B^ — > and B$ 3> A^ is identical but with $ and \I/ interchanged in the 
expressions above. In both cases, the steady-state particle populations are on the order of 
NiB t § ~ 10 25 and ~ 10 31 (when they are non-zero) for the fiducial input values listed 

above. 

The time needed to reach this steady state from a negligible initial dark matter density 
in the star is on the order of 

t ss ~ max j 5 " 1 ' • ( 31 ) 

where the first value corresponds to IND balancing capture and the second to annihilation 
balancing production of anti-hDM by IND. For the fiducial cross sections we are considering, 
we find that annihilation takes longer to balance, and leads to t ss ~ 2 x 10 7 s. This is ultra 
short relative to the lifetime of a typical neutron star. 

Having reached a steady state, the main combined effect of the IND and annihilation 
processes is to inject energy into the host neutron star with rate (mj + m$ + mjv)C. In 



For Ai > Bi, we find: Nq, = AT$ = C/(B* + B$), N$ = B$/A 9 , = B^/A^ 
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Figure 6: The build-up of the abundance of hDM and anti-hDM in a neutron star for 
mq, = 2.85 GeV and m$ = 2.05 GeV for a total dark matter density of 0.3 GeV/cm 3 . At 
this density a steady state is attained after a little over 10 7 s. 

this respect, hDM in the steady state limit has the same effect on neutron stars as ordinary 
self-annihilating DM. Energy injection by DM can interfere with and halt the cooling of old 
neutron stars, and therefore the observation of a very cool, old neutron star in a region of large 
DM density could put significant constraints on a wide variety of DM scenarios [HJ SHI HH] ■ 
However, this effect is too small to be observed using existing observations, and appears 
to be challenging to probe in the near future [HJ HH1 US]. With hDM, baryons within the 
neutron star are also destroyed by IND, but the number is negligible compared to the total 
of Nb — 2 x 10 57 within a typical neutron star over the lifetime of the Universe unless the 
local DM density approaches an enormous value of 10 14 GeV/cm 3 . 

Consider next the case of a vanishingly small IND cross section. The populations of 
of \I/ and $ hDM particles will now build up within the neutron star with rate given by 
Eq. ( 126|) . corresponding to populations of about lO* 3 (pdm/ GeV cm -3 ) over the lifetime of 
the Universe. We can compare this number to the populations required for DM particles to 
begin self-gravitating and to form a black hole. Self-gravitation begins when [HJ |50l [51] 

1.4 x 10 18 kg/m 3 V 2 / „„, 



This is much greater than the steady-state populations found above for the case of a large 
IND rate, and also greater than the numbers collected for a smaller IND rates for local hDM 
densities below (3 x 10 2 GeV/cm 3 ) min{l, 3 x 10 _57 cm 3 s _1 /( crt ')-fA r £ ) }- 

Self-gravitating DM particles will form a black hole that can potentially destroy the host 
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star unless they are stabilized in some way jHJ H7J EH EU [52] . For non- interacting fermions 
there is a degeneracy pressure that must be overcome. To do so, the total number of fermions 
must exceed 

3 



3 

(33) 



with Mpi = y^8n/G ~ 2.4 x 10 18 GeV the reduced Planck mass. In the case of non-self- 
interacting bosons, there is still a zero-point pressure which can be overcome if the number 
of bosons exceeds 

2 

- ' ] ~2x 10 37 

rrii I \ rrii 



N.ZN^l^l)'^,^^) 2 . ,34) 



A bosonic black hole can also arise even before bulk self-gravitation sets in through the 
formation of a Bose-Einstein condensate 



We argue that black holes are not likely to form in the specific theory of hylogenesis 
discussed in Ref. [To] and Appendix |AJ in which $ and \1/ carry equal and opposite charges 
under a spontaneously broken U{1)' gauge symmetry, until Ny t $ > iV"/ rii . At distances 
much smaller than the inverse vector mass m~ z } ) the effects of breaking the U{1)' can be 
neglected and the Z' vector boson mediates a repulsive force between particles with like- 
sign charges. This induces a pressure among $ particles that prevents them from collapsing 
into a black hole provided mz' <C (m$M|j) 1//3 and e' ^> m$/Mp\ (where e' is the U(l)' 
coupling of $) [50J. Instead, under these mild assumptions the formation of a black hole 
requires an approximately charge-neutral collection of hDM particles and therefore roughly 
equal numbers of \l/ and $ states. A necessary condition for this to occur is N > N^ rit to 
overcome the fermion degeneracy pressure. Therefore we do not expect this specific theory 
of hylogenesis to lead to the formation of hDM black holes within neutron stars unless the 
local hDM density around the host star exceeds 5 x 10 11 GeV/cm 3 . 

Between the two extremes of the large and small IND rates considered above, no new 
obvious observational bounds arise. With a small but non-zero rate for IND, the \l/ and $ 
populations may grow large enough to form some amount of anti-hDM that will subsequently 
annihilate away. Large hDM densities can also lead to direct annihilation of hDM to 
antibaryons, ^ $ — > N M, and similarly for anti-hDM. This can be accounted for by adding 
terms of the form —DN^N^> to the evolution equations for \l/ and $, with D ~ (cv)iND/Vf l DM 
and VhDM is the volume occupied by the hDM, whose evolution we must also keep track of. 
Note that VhDM appears here rather than 47irf h /3 because this effect starts to compete 
with standard IND only after the hDM population begins to self-gravitate, when iVj > 
PcV}iDm/V > N se if. We defer an analysis of the effects of self-gravitation and self-annihilation 
of hDM on the structure of neutron stars to a future work. 



4.3 White Dwarfs 

We can perform a similar analysis for the effects of hDM on white dwarfs. A typical white 
dwarf consists primarily of carbon and oxygen and is supported by the degeneracy pressure 
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of the electrons it contains. Typical white dwarf parameters are R = 0.01 R G} M = 0.7 M , 
p c = 10 9 kg/m 3 , T c = 10 7 K [53j. We will approximate the internal structure of a white dwarf 
as consisting entirely of carbon with a uniform density 



The capture rate of DM through elastic scattering with nuclei in a white dwarf is 
approximated well by the expression [UJJ HU [5U [55j EH] 

with v being the local velocity dispersion, v esc the escape velocity at the surface of the star, 
and the the effective cross section is defined to be 



a e ff = mm 



Ji 1" ) I (36) 



k 




where the sum runs over all nuclear species k in the star with mass fraction Xk and atomic 
mass (number) (Zk), a^ 1 is the spin- independent nucleon scattering cross section off 
protons, f p and f n are the relative coupling strengths of the DM to protons and neutrons, 
and m Tk {m rp ) is the reduced mass of the DM- nucleus (DM-proton) system. The saturation 
of the effective cross section at ttR 2 corresponds to the star becoming optically thick. 

We find that a typical white dwarf is indeed optically thick for the fiducial scattering 
cross section of a^ 1 = 10~ 39 cm 2 (with f p = 1, f n = 0). This gives 

= C$ = C (37) 

- (6 x lO^s- 1 ) (— M f — ( PhDM \ ( 5GeV ) f 270km /s \ 

where R & ~ 7.0 x 10 5 km is the solar radius. Once captured, hDM particles thermalize 
rapidly and collect within a thermal radius at the core of the star of size 



rr in y s ^3 GeV\ l/2 / T c \ 1/2 / 10 9 kg/m sX ' " 



(38) 



For my = 771® = 2.25 GeV, we find A { = 1 x lO" 4 ^" 1 and B i = 6x lO^V 1 (Os^ 1 ) for the 
case of a large (small) IND rate (<jv)ind = 10~ 39 cm 3 /s (0cm 3 /s) 

In the case of a large IND rate, we find that the populations of $, *Sf, and $^ all reach 
a steady state quickly relative to the lifetime of a typical white dwarf. In this state we find 
iV,j, ~ ~ 5 x 10 36 and iVq> ~ ~ 6 x 10 49 for the fiducial parameter values listed above. 
The number of baryons destroyed by IND over the age of the Universe is a small fraction 
of the total number in a white dwarf provided pdm "C 10 11 GeV/cm 3 . The main effect of 
hDM capture, IND, and annihilation on white dwarf is therefore the injection of energy into 
the stellar interior with rate (m^, + m$ + mjv)C. This is equivalent to the effect of ordinary 
self-annihilating DM (with slightly larger mass). 



Heating of white dwarfs by DM capture and annihilation was studied recently in Refs. 
[56] . These authors differ in their conclusions, with Ref. [55] finding an upper bound on 
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the nuclear capture cross section of about 10~ 43 cm 2 and Ref. [56] finding essentially no 
bounds from current observations. The origin of the disagreement is the density of DM 
within globular clusters, from which Ref. [55] obtains their most stringent bounds, while 
Ref. [56] argues that they contain a much smaller dark matter abundance. The observation 
of a cool white dwarf within a dwarf spheroidal galaxy could help to resolve this question 
decisively [56]. Based on this uncertainty, we do not consider hDM models with nuclear 
scattering cross sections as large as cxf 7 = 10~ 39 cm 2 to be ruled out. However, we do note 
that the model of hDM considered here can also accommodate much smaller elastic scattering 
cross sections. 

With a small IND rate, hDM collects and builds up within the stellar interior. As 
we argued above, we do not expect any significant effects on the star until at least Ni ~ 
Nseifi ^Lit gi ven m Eqs. (I32|33p . For this to occur over the lifetime of the Universe, a local 
hDM density approaching 10 8 — 10 10 GeV/cm 3 is needed. We do not know of any observations 
of white dwarfs that approach these requirements. 



4.4 The Sun and other Main- Sequence Stars 

The capture rate of light (1—5 GeV) DM in the sun is dominated by scattering with hydrogen 
and helium and is given to a good approximation by [101 [57] 

c . „ (g x lo2V1) ( 5GeV \ / PDM \ / 270km/s \ / _of 
1 ~ 1 '{my + mzj V0.3GeV/cm 3 y' V v )[ 10" 39 c 



cm 2 



2 

TfL 

x H + (l.l)x He (l + f n /f p ) 2 ^ 

mi 

1 v 



(39) 



where the relative factor for helium relative to hydrogen comes from its slightly different 
distribution within the sun. In writing this expression we have also made use of the fact that 
the sun is optically thin for our fiducial nucleon elastic scattering cross section. A second 
consequence of the sun being optically thin is that nearly all the anti-hDM produced by IND 
escapes from the sun, since it is typically produced with velocities much larger than v esc . 
Thus we set $t — >• in our evolution equations. The thermal radius in which the DM 
collects is 

Tim . (5 x !0» c m > " 2 ( U> "V-' ) " 2 ,40) 



mi J \,1.5xl0 7 K/ V p, 

which yields Bi ~ 10 _13 s _1 (0s _1 ) for the large (small) IND rate. 

An additional effect that is important for the sun is evaporation, in which captured hDM 
particles are up-scattered above the escape velocity and leave the sun [391 US]- Fitting to 
the results of Refs. [15J |57] , we take the evaporation rate to be 

E- ~ in[-3-6(WGIeV)-4] [ ^ \ -1 (A] \ 

l ~ U 5xlO- 3 W ' 1 ' 
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This evaporation adds a dissipative —EiNi term to each of Eqs. (!20] - [23l) . For the fiducial 
parameters we are considering (and a large IND rate), evaporation becomes more important 
than IND for hDM masses below rrii < 2.4 GeV. 

We find that capture, evaporation, and IND reach a steady state on a timescale of 1CT 5 - 
lCT 1 billion years, much shorter than the age of the sun. The equilibrium populations are 
always less than iVj ~ Ci/max{B il E{\ < 10 41 , for both large and small IND rates. This is a 
very small fraction of the total mass of the sun, and is far below what is needed for hDM to 
self-gravitate or to contribute significantly to the properties of the solar interior [58] . Note 
as well that due to the low masses of hDM, the energies of any neutrinos produced by IND 
are well below the thresholds of DM searches for upward-going muons in neutrino telescopes 
such as Ice Cube [571 [59]. The effects of hDM on other main-sequence stars are also expected 
to be small, unless the star is immersed in a region of very large DM density [60] . 

5 Conclusions 

We have investigated signals of nucleon destruction that can arise in specific theories where 
dark matter is antibaryonic. Antibaryonic DM is motivated by unified mechanisms for DM 
and baryon generation (referred to as hylogenesis) which address the cosmic coincidence 
between the energy densities of dark and visible matter in the Universe. 

In the hylogenesis scenario considered here, the DM consists of an asymmetric density 
of fermions ^ and scalars $. These particles can scatter inelastically with nucleons via 
the reactions ^ N — > $t M and $ iV — > M, where N is a nucleon and M is a meson. 
These induced nucleon decay (IND) processes lead to distinctive signatures in nucleon decay 
searches, stellar evolution, and hadron colliders. 

IND is a novel signature of DM that can be searched for in terrestrial nucleon decay 
searches. The effective nucleon lifetime is expected to be ~ 10 29 — 10 32 years, if baryon 
transfer between dark and visible sectors is mediated by new physics at a scale A ~ 300 GeV— 
1 TeV. Since the DM states are unobserved, IND events mimic standard nucleon decay 
with neutrino final states, but typically with greater final state meson energy E ~ GeV. 
Due to these different kinematics, existing searches do not directly apply in general. Our 
study therefore motivates new searches in these experiments. We expect that the resulting 
sensitivities should be comparable to those in standard nucleon decay searches. 

The coupling of hylogenic DM to quarks responsible for IND can also give rise to 
observable signals in hadron colliders. In particular, such couplings can potentially lead 
to monojet signals at the Tevatron and the LHC. Existing monojet searches at the Tevatron 
place a lower bound on the coupling of hylogenic DM to quarks, and this bound will be 
significantly improved at the LHC. We find that the coupling strengths that can be probed 
at the LHC are of the same size as those that can produce an observable effect in nucleon 
decay search experiments. This correlation may permit the characterization of hylogenic 
DM through a very diverse set of experimental probes. 

If hylogenic DM also scatters elastically with nucleons, it can be captured in stars. Once 
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captured, it will thermalize with the baryons in the star and sink to the stellar core. In 
compact stars, the local baryon density in the stellar core may be large enough for IND to 
occur at a significant rate, destroying both a nuclon and the hDM particle, and producing 
an anti-hDM particle and a meson. This process deposits energy in the star, an effect that 
may be probed by observing old and cool neutron stars or white dwarfs. However, current 
observations do not provide a definitive bound on or evidence for these processes. 
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A Appendix 

In this Appendix we review some of the important features of the hylogenesis model presented 
in Ref. [15J. The model consists of the SM together with a SM-neutral hidden sector 
containing two massive Dirac fermions X a (a = 1,2, with masses rrix 2 > rnx 1 > TeV), 
a Dirac fermion and a complex scalar $ (with masses m# ~ m$ ~ GeV). These fields 
couple through the "neutron portal" (XU C D C D C ) and a Yukawa interaction: 

- Sf D |p «i4)K4) + Ca (X a , L * L + X a , A)$ + h.c. (42) 

where i,j, k label flavors and the quark color indices are implicitly contracted antisymmet- 
rically. Many variations on these operators exist, corresponding to different combinations of 
quark flavors and spinor contractions. With this set of interactions one can define a conserved 
generalized global baryon number with charges Bx = —(B^, + -B$) = 1. The proton, and 
$ are stable due to their B and gauge charges if their masses satisfy 

— m$| < m p + m e , m p — m e < + m$ . (43) 

\l/ and $ are the "hidden antibaryons" that comprise the dark matter. Furthermore, there 
exists a physical CP-violating phase arg(A*A2CiC:2*) that cannot be removed through phase 
redefinitions of the fields. 



21 



We also introduce a hidden U(l)' gauge symmetry under which \1/ and $ have opposite 
charges ±e', while X a is neutral. We assume this symmetry is spontaneously broken at the 
GeV scale, and has a kinetic mixing with SM hypercharge U (l)y via the coupling —^B^Z', 
where and Z' are the U(l)y and U(l)' field strength tensors. At energies well below 
the electroweak scale the effect of this mixing is primarily to generate a vector coupling of 
the massive Z' gauge boson to SM particles with strength — cwKQ em e. The GeV-scale Z' 
masses we consider here can be consistent with observations for 10~ 6 < k < 1CT 2 [61, 62J. 

In the mechanism for hylogenesis presented in Ref. [15], baryogenesis begins when a 
non-thermal, CP-symmetric population of X% and X\ is produced in the early Universe. 
These states decay through X\ — > udd or X\ — > (and their conjugates). An asymmetry 
between the partial widths for X\ — > udd and X\ — > udd arises from interference between 
tree and loop diagrams and is characterized by 

e = — [r(A x udd) - T( Xl udd)} ~ ^^iy 2 ^ 21 , (44) 
2T Xl 1 v ' v Ji 256vr 3 |Ci | 2 M 4 m X2 ' v ; 

where we have assumed that the total decay rate Tx t is dominated by X\ — > over the 
three-quark mode, and that mx 2 ^> rnx 1 - For e ^ 0, X± decays generate a baryon asymmetry 
in the visible sector, and by CPT an equal and opposite baryon asymmetry in the hidden 
sector. These asymmetries can be "frozen in" by the weakness of the coupling between both 
sectors provided the temperature at which the X x are produced is not too high. For the 
asymmetry to be large enough to explain the observed value, mx 1 2 cannot be too much 
smaller than M. 

Once produced, \& and $ will thermalize by scattering with Z' vectors present in the 
plasma. These interactions will also deplete the symmetric densities of ^ and $ very 
efficiently through annihilation to pairs of Z' vectors provided mz> < [HS]- Only the 

asymmetries will remain. This is analogous to the annihilation of baryons with antibaryons. 
The cross section for tyty — > Z'Z' is given by 



Annihilation of is given by a similar expression. 

There is also a direct detection signal in our model due to the hidden Z' mediating the 
elastic scattering of ^ and $ off protons. The effective scattering cross section per nucleon 
for either ^ or $ is spin-independent and given by 

where is the DM-nucleon reduced mass. For a DM mass of 2.9 GeV, this is slightly below 
the best current limits from CRESST @T], CDMS @2], and CoGeNT g3]. 
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